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ABSTRACT 

Damped Lyman-alpha (DLA) and sub-DLA absorbers in quasar spectra provide the 
most sensitive tools for measuring element abundances of distant galaxies. Estimation 
of abundances from absorption lines depends sensitively on the accuracy of the atomic 
data used. We have started a project to produce new atomic spectroscopic parame¬ 
ters for optical/UV spectral lines using state-of-the-art computer codes employing very 
broad configuration interaction basis. Here we report our results for Zn II, an ion used 
widely in studies of the interstellar medium (ISM) as well as DLA/sub-DLAs. We 
report new calculations of many energy levels of Zn II, and the line strengths of the 
resulting radiative transitions. Our calculations use the configuration interaction ap¬ 
proach within a numerical Hartree-Fock framework. We use both non-relativistic and 
quasi-relativistic one-electron radial orbitals. We have incorporated the results of these 
atomic calculations into the plasma simulation code Cloudy, and applied them to a lab 
plasma and examples of a DLA and a sub-DLA. Our values of the Zn II AA 2026, 2062 
oscillator strengths are higher than previous values by 0.10 dex. Cloudy calculations for 
representative absorbers with the revised Zn atomic data imply ionization corrections 
lower than calculated before by 0.05 dex. The new results imply Zn metallicities should 
be lower by 0.1 dex for DLAs and by 0.13-0.15 dex for sub-DLAs than in past studies. 

Our results can be applied to other studies of Zn II in the Galactic and extragalactic 
ISM. 

Subject headings: atomic data; atomic processes; ISM; abundances; Galaxies: abun¬ 
dances; quasars: absorption lines 


1. Introduction 

Most elements heavier than He are produced by stellar evolution and then distributed into 
interstellar space. Understanding the chemical composition of distant galaxies is therefore crucial 
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to understanding the star formation and feedback processes central to galaxy evolution. Absorption 
lines of DLA and sub-DLA absorbers in the spectra of background quasars provide the most sensitive 
tools to measure the heavy element content of distant galaxies. The DLAs have neutral hydrogen 
column densities A^hi > 2 x cm“^, and the sub-DLAs have 10^® < A^hi < 2 x cm“^. The 
DLAs and sub-DLAs dominate the neutral gas content of galaxies, and constitute most of the H I 


available for star formation at redshifts 0 < z < 5 [e.g., Storrie-Lombardi &: Wolfe (2000); Peroux et 


al. (2005); Prochaska et al. (2005); Noterdaeme et al. (2012); Zafar et al. (2013)]. DLAs observed 


toward GRB afterglows also offer an excellent probe of the physical and chemical conditions in 
distant galaxies [e.g., 


Savaglio et al. 

(2003) 

Chen et al. 

(2005) 

Prochaska et al. 

(2007) 

Fynbo et 


al. (2009)]. 


Besides neutral hydrogen, DLAs and sub-DLAs also show a number of other elements ranging 
from C to Zn. The abundances of these elements provide very sensitive indicators of the chemical 


evolution of galaxies [e.g., 

Pettini et al. 

(1997 

); 

Kulkarni Sz Fall 

(2002 

); 

Prochaska et al. 

(2003 

); 

Peroux et al. 

(2008) 

Meiring et al. 

(2009b) 

Cooke et al. 

(2011 

); Rafelski et al. 

(2012) 

Som et al. 


(2013, 2014)]. Since different elements are produced by stars of different masses, the measurements 
of element abundances as a function of time give information about the history of formation of 
stars of different masses in galaxies. 

The quality of the atomic data directly affect the accuracy of the element abundances and 
physical properties of galaxies that are estimated from the measurements of absorption lines. The 
most commonly used atomic data reference for the analysis of absorption lines in DLAs and sub- 


DLAs is Morton ( 

2003) [see, e.g., Battisti et al. 

(2012); Rafelski et al. 

(2012); Kulkarni et al. 

(2012); Cuimaraes 

(2012); Jorgenson et al. (2013) 

Som et al. 

2013 2014 

)]. On the other hand. 

there is a need to improve beyond the oscillator strengths of 

Morton (2003). This is because 


Morton (2003) lists large uncertainties for the oscillator strengths of some important transitions. 
[These uncertainty values are listed on the NIST Atomic Spectra Database ( Kramida et al.|[2014 )]. 
Furthermore, for some transitions, Morton (2003) lists no oscillator strengths at all. In some 


cases, the NIST database assigns low accuracy grades even for more recent values obtained since 


Morton (2003). Such shortcomings in atomic data limit our ability to interpret the spectra of 


high-redshift galaxies, potentially leading to erroneous inferences about their chemical enrichment 
and star formation history. 


With the goal of producing new reliable atomic data for commonly used astrophysical ions, we 
have started a collaborative study combining atomic physics, plasma simulations, and observational 
spectroscopy. The goals of this study are to examine the quality of available atomic data, to improve 
the accuracy of the atomic data with low reported accuracies, to incorporate these new atomic data 
into Cloudy (our widely utilized plasma simulation code), and to study the effect of the revised 
Cloudy code on the analysis of absorption lines in DLAs and sub-DLAs. In a recent study (Kisielius 
et al. 2014), we examined the atomic data for the key ion S II. Here we focus on the ion Zn II, 
which also plays a very important role in studies of DLAs and sub-DLAs. 
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1.1. Why Zn II? 


Refractory elements such as Fe, Si, or Mg condense in the form of solid dust grains in the 
interstellar medium. On the other hand, volatile elements such as N, O, P, S, Ar, Zn do not 
condense appreciably on interstellar dust grains. The gas-phase abundances of such weakly depleted 
elements can therefore give their total (gas -|- solid phase) abundances. For observations of DLAs 
and sub-DLAs, weak, unsaturated lines of the elements N, O, P, or Ar are often not accessible 
in ground-based spectroscopy, which makes it difficult to measure their column densities reliably. 
Moreover, the abundance of N appears to be complicated by nucleosynthetic differences between 
primary and secondary N production [e.g., Pettini et al. (1995); Som et al. (2014)]. The two 
elements that have emerged as the most useful for DLA/sub-DLA studies are S and Zn. Having 
discussed the atomic data for S in Kisielius et al. (|2014), we now turn to Zn. 


Zn is especially interesting because it tracks Fe closely in Galactic stars (for [Fe/H] > —2). 
Especially important among the Zn ions is Zn II, which is the dominant ion in DLAs. A key 
advantage of Zn II is that it has two weak absorption lines in a relatively narrow wavelength re¬ 
gion at A A2026.14, 2062.66, that are often unsaturated and hence allow accurate column density 
determinations. Moreover, these lines often lie outside the Lyman-a forest, and therefore allow un¬ 
ambiguous identifications and measurements free of blends. This has made Zn the most commonly 


used metallicity indicator for DLAs and sub-DLAs. Starting from the early work of Meyer (1989) 


and Pettini et al. (1994), Zn has been observed in > 150 DLAs and ~ 50 sub-DLAs at redshifts 


ranging from z < 0.1 to z > 3.3. Zn has also been used as a metallicity indicator in studies of the 
Milky Way interstellar gas. For all of these reasons, accurate atomic data for Zn II transitions are 
very important. 


Morton (2003) lists the oscillator strengths of 0.501 and 0.246 for the Zn II A A 2026.14, 2062.66 


lines, respectively, but there are no estimates of the uncertainties in these values in the NIST 
database. Thus, the uncertainties in the metallicity introduced by the uncertainty in the oscillator 
strengths could be far larger than those often quoted from the measurement uncertainties in high- 
resolution data (typically < 0.05 dex). Zn II has additional absorption lines at 923.98, 938.71, 
949.46, 984.14, and 986.52 A, but they are listed in Morton (2003) without any oscillator strength 
estimates. 


1.2. Previous Zn II Calculations 


There is quite a generous amount of either experimental or theoretical studies considering the 
transition wavelengths, radiative transitions or scattering processes in the ion Zn II, see, e.g., exper- 


imental works of 

Bergeson & Lawler 

(1993); 

Mayo et al. 

(2006) 

; Gullberg &: Litzen 

(2000 

). Multi- 

configuration Hartree-Fock calculations were reported by 

Froese Fischer 

(1987 

). The Hartree-Fock 

approximation adopting transformed radial orbitals was used by 

Karpuskiene & Bogdanovich 

(2001) 


to calculate oscillator strengths of astrophysically important lines in Zn I and Zn II ions. Their cal- 
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culations which included core-polarization effects produced oscillator strengths which agreed quite 


well with previously published semiempirical values. Recently Harrison & Hibbert (2003) have 


presented a pivotal study of oscillator strengths for the Zn II. They have investigated the 4s-4p 
resonance line oscillator strengths in the Zn II ion using extensive configuration interaction (Cl) cal¬ 
culations. They studied the influence of core-polarization, electron-correlation in the core, core-core 
correlation effects in resulting oscillator strengths using CIVS computer code which deals with rela¬ 
tivistic effects in Breit-Pauli approximation. Their determined oscillator strength values, / = 0.268 
for the 4s ^S'i/2 — / = 0.547 for the 4s^S'i/2 —line he about 5 — 10% higher 

than the recent experimental values obtained by Bergeson fc Lawler] (1993), / = 0.255 ± 0.024 and 


/ = 0.492 ± .039, respectively. Nevertheless, the theoretical / values of Harrison &: Hibbert (2003) 


are in good agreement with relativistic many-body perturbation theory calculation results of Chou 


&: Johnson (1997), which are / = 0.264 and / = 0.538 for these lines. Unfortunately, Harrison 


&: Hibbert! (2003) have considered only 4s - 4p lines, therefore their data are not enough for a 


comprehensive modeling calculations and can not be employed in our investigation. 


Very recently (()elik et ah (2013) published new calculations of atomic data for the Zn II ion. 


Their calculations were performed using two different semi-empirical methods, the weakest bound 
electron potential model theory (WBEPMT)and the quantum defect orbital theory (QDOT). They 
employed numerical Coulomb approximation wave functions and numerical non-relativistic Hartree- 
Fock wave functions to determine the necessary parameters. As a result, the real multi-particle 
system was transformed to a simple one-particle system. In the WBEPMT case, the effective 
parameters (nuclear charge Z*, n*, I*) were derived from the experimental energy data or from 
other existing calculations for the mean radius < r >. One can have a reasonable doubt if these 
< r > values, derived from the different sources, are really coherent and consistent. Unfortunately, 
their suitability for the 4s electrons is very questionable as their wavefunctions significantly overlap 
with those of other electrons. As a result, the results from these two different approximations differ 
quite significantly one from another even if they both are semi-empirical and based on the same 
experimental data. 


Several works have considered electron-impact excitation parameters for the Zn II ion. Scatter¬ 


ing parameters for some transitions in this ion were determined by Pindzola et al. (1991); Zatsarinny 


&: Bandurina (1999) by employing R-matrix methods in the LS'-coupling approximation. Sharma 


et al. (2011) have applied a fully relativistic distorted-wave theory to study the electron-impact ex¬ 


citation of the ns — np resonance transitions in singly-charged metal ions with one valence electron, 
including Zn"*" ions. Unfortunately, the above-mentioned studies did not produce atomic data sets 
suitable for spectral modeling where the complete and consistent data are required. 


In order to assess the accuracy of the Zn H atomic data, we performed new calculations of the 
oscillator strengths for all Zn II electric dipole, magnetic dipole, and electric quadrupole transitions. 
Section 2 describes our new calculations and how they compare with previous estimates. Section 
3 outlines the inclusion of these calculations into Cloudy and gives a few examples of applications. 
Finally, section 4 summarizes our results and their implications for abundance studies of DLAs and 
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sub-DLAs. 


2. Calculations of new atomic data 


Our calculations are performed by employing Hartree-Fock radial orbitals (HFRO). The rela¬ 
tivistic corrections are inclnded in the Breit-Pauli approximation. We determine spectral parame¬ 
ters for four even configurations 3d®4s^, ,and 3d^^5s and for three odd configurations 

Sd^^Ap, 3d^^5p, and Sd^AsAp. The configuration Sd^AsAp levels lie in an energy range which is signif¬ 
icantly wider compared to the purpose of this work. For that reason we determine only the lowest 
levels of this configuration arising from the 4s and Ap electrons bound into term. The electron- 
correlation effects are included in the configuration interaction (Cl) approximation by adopting the 


basis of transformed radial orbitals (TRO) described by Bogdanovich &: Karpuskiene (1999). 


At the first step, we solve the HF equations for the ground configuration using the code of 


Froese Fischer (1987). Next, we determine solutions of HF equations P{nl\r) for all Al, 5s and 5p 
orbitals in a frozen-core potential. This basis of HFRO is complemented with transformed radial 
orbitals PTRo{nl\r), which are introduced to describe virtual electron excitations from adjusted 
(investigated) configurations. The TRO are obtained by a way of transformation: 


PTRo{nl\r) =N{P '“+''exp(-Hr)P(no/ok) 

poo 

— P{n'l\r) / P{n'l\r')r'exp{—Br')P{nolo\r')dP) 


( 1 ) 


n'<n 


The parameters k and B are introduced to ensure the maximum of the mean energy corrections 


determined in the second-order perturbation theory (PT) (see Bogdanovich &: Karpuskiene (1999)). 
Here the factor N ensures the normalization of the TROs, which are determined for the electrons 
with the principal quantnm number values 6 < n < 11 and for all allowed values of the orbital 
quantum number 1. 

The set of admixed configurations is generated from the adjusted configurations by introducing 
one-electron and two-electron excitations from the 3p, 3d, 4s, 4p, 4d, and 5s shells to all available 
states in the basis of determined radial orbitals. This leads to a huge set of admixed configurations 
and, consequently, to large Hamiltonian matrices. In order to reduce the size of the Hamiltonian 
matrices to be diagonalized, we need to determine the most significant configurations. As a se¬ 
lection criterion, we use the averaged weights W of the admixed configurations K'T'LS in the Cl 
wavefnnction expansion of the adjusted configuration KqTLS: 


Wrt{K K') + 1)(2S + l){KoTLS\\H\\K'rLS)^ 

g{Ko) {E{K') - E{Ko)f 


( 2 ) 
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where T describes all possible intermediate momenta, which bound the non-relativistic configura¬ 
tions K and Kq into LS multiplet. These averaged weights are determined in the second-order PT. 
We select only those configurations which have their weights larger than some selection criterion 
Wq. More details on configuration selection procedure are given by Bogdanovich &: Karpuskiene 

pooTl). 


Bearing in mind restrictions of our computer system, we can perform our spectroscopic data 
generation for Wq ~ 5 x 10“®. The main limiting factor is the size M of the electrostatic interaction 
operator Hq which depends on the number of configuration state functions (CSF) having the same 
total momenta LS. We have performed two other similar test calculations with the value of Wq 
increased by a factor of two (CIJ^^q) and by a factor of ten (CI!J|[q). In the further description 
we present the parameters of these three calculations separated by a slash. In that way we can 
demonstrate the convergence of our calculated results. 


For the four even configurations investigated, we have selected the most important = 
1012/745/338 non-relativistic configurations, including adjusted ones, which give rise to Ce = 
749470/517350/209921 configuration state functions (CSF). Here we apply CSF reduction by mov¬ 
ing the shells of virtually excited electrons to the beginning of active shells, as it was described 
in detail by Bogdanovich et al. (2002); Bogdanovich &: Momkauskaite (2004). This procedure re¬ 
duces the number of CSFs to Re = 20991/16021/7522. Using this basis, the size of the largest 
Hq matrix is Me = 17317/13195/6131. For the three odd configurations investigated, we have 
selected So = 784/577/252 most important configurations, including the adjusted ones. These con¬ 
figurations make Co = 2513717/1727977/721017 CSFs, which in turn are further reduced by the 
CSF-reduction procedure to Ro = 531815/386459/181535. Adopting this base, the size of largest 
Hq matrix reaches Mq = 108720/78903/36696. 


The Hamiltonian eigenvalues and eigenfunctions are determined adopting this reduced configu¬ 
ration basis. Further, the determined Cl wavefunctions are employed to calculate electron transition 
parameters. The Ml and E2 radiative transition data are produced for transitions among the levels 
of the same-parity configurations, and the El, M2, and E3 transition parameters are produced for 
the radiative transitions among the levels of the different-parity configurations. The significance of 
the radiative transitions of higher multipole order, such as M2 or E3, for the radiative lifetimes of 



tation collision strengths in the plane-wave Born approximation. The most inclusive description 


of the adopted approach is given by Bogdanovich 

(2004 

2005 

), whereas its application for data 

production and their accuracy analysis is given by ’ 

Bogdanovich Sz Martinson ( 

1999) 

Bogdanovich 


et al. 

(2003a|b); 

Karpuskiene &: Bogdanovich 

(2003) 

Kupliauskiene et al. 

(2006 


To perform our calculations, we have employed our own computer codes together with the 
codes of Hibbert et al. (1991); Froese Fischer et al. (1991); Froese Fischer &: Godefroid (1991) 
which have been adopted for our computing needs. The code of Hibbert et al. (1991) has been 
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updated according to the methods presented by Gaigalas et al. (1997 1998). 


2.1. Energy levels and wavelengths 


We present a comparison of our calculated Zn 11 level energies with the data extracted from the 
NIST database (Kramida et al. 2014) in Table As we have explained in Sect. three different 
Cl expansions are adopted in the Zn II level energies calculation using the same multiconfiguration 
Hartree-Fock method based on the TRO. As one can clearly see, the agreement of the data becomes 
closer when the selection criterion Wpt goes down. The mean-square deviations MSD are given 
in the last row of Table for the indication of the convergence of our calculation. One can notice a 
considerable improvement of accuracies of the calculated energies when the selection criterion Wpt 
is decreased two times (from MSD = 3809 to MSD = 1146 cm“^), and a smaller improvement 
when this criterion is further decreased five times (down to MSD = 730 cm“^). 

Alongside the level energies determined using the different Cl expansions, we present the 
corresponding percentage deviations 6E: 


5E = - Enist . (3, 

F/NIST 

As one can notice, the 6E decrease consistently for most levels, as the Cl expansion increases. This 
again is underlining the fact that (i) the accuracy of our energy levels increases when the extended 
Cl basis is employed; and (ii) the convergence for the calculated level energies is achieved as the 
final changes in the 5E values are much smaller compared to the initial ones. Based on our past 
experience with such calculations, we do not expect further extensions of the Cl basis to yield any 
substantial changes of the level energies or the radiative transition parameters. 


Here we should explain that in further consideration or in modeling performed by Cloudy we 
do not use the calculated energies. They are substituted by the experimental values in order to 
enable us to have the correct transition wavelengths. The main reason to present Table here 
is to demonstrate the convergence of our calculations for the energy levels, which reflects on the 
production of transition parameters such as oscillator strengths gf or transition probabilities A. 


Using the determined Cl wavefunction expansions, we are able to generate data sets for the 
radiative transitions involving twenty seven energy levels. In Fig. we compare the weighted 
oscillator strengths determined in approximation CItro with gf values determined in CF^^q and 
CF^^q approximations. Although there are some weaker lines where deviations can exceed 25%, 
especially in CI^|[q approximation, for most lines the agreement is within the range of 25%. One 
can notice that the deviations have become significantly smaller in CF^j^q approximation, indicating 
the convergence of the gf values, as it has been a case with level energy values (see Table [^. 

Table gives a sample of the transition parameters used in the plasma simulation package 
Cloudy. Here transition line strengths S are tabulated together with the type of transition (El, 











Table 1. 


Comparison of calculated Zn II level energies E (in cm and their percentage 
deviations 5E with experimental data from the NIST database. 


N 

State 

2J + 1 

NIST 

CItro 

dPxRO 

i^jredl 

'^AxrO 

c Tjiredl 

^-t^TRO 

^Tred2 

'^J^TRO 

r 77'red2 

^-t^TRO 

1 

3di°4s 

2 

0.00 

0 

- 

0 


0 

- 

2 

3d^°4p 

2 

48481.00 

48670 

0.39 

47530 

-1.96 

48546 

0.13 

3 

3d^°4p 

4 

49355.04 

49144 

-0.43 

48004 

-2.74 

49018 

-0.68 

4 

3d®4s2 2p 

6 

62722.45 

62815 

0.15 

63515 

1.26 

66605 

6.19 

5 

3d®4^2 

4 

65441.64 

65404 

-0.16 

66110 

1.02 

69292 

5.88 

6 

3d^°5s 

2 

88437.15 

88394 

-0.05 

88059 

-0.43 

92729 

4.85 

7 

3d^°4d 

4 

96909.74 

96972 

-0.06 

96864 

-0.05 

102361 

5.62 

8 

3d^°4d 

6 

96960.40 

96988 

0.03 

96880 

-0.08 

102377 

5.59 

9 

3d^°5p ®P 

2 

101365.9 

101571 

0.20 

100818 

-0.54 

106727 

5.29 

10 

3d^°5p ®P 

4 

101611.4 

101705 

0.09 

100953 

-0.65 

106865 

5.17 

11 

3d^(^D)4s4p*P 

6 

103701.6 

105603 

1.83 

106286 

2.49 

108677 

4.80 

12 

3d^l^D)4s4p*P 

4 

105322.7 

106855 

1.45 

107526 

2.09 

109935 

4.38 

13 

3d^l^D)4s4p*P 

2 

106528.8 

107849 

1.24 

108530 

1.88 

110936 

4.14 

14 

3d®(®P)4s4p^P 

10 

106779.9 

107465 

0.64 

108014 

1.16 

110474 

3.46 

15 

3d^(^D)4s4p*F 

8 

106852.4 

107938 

1.03 

108470 

1.51 

110940 

3.82 

16 

3d^CD)4s4p*F 

6 

107268.6 

108515 

1.16 

109039 

1.65 

111522 

3.96 

17 

3d^l^D)4s4p*F 

4 

108227.9 

109435 

1.11 

109973 

1.61 

112437 

3.89 

18 

3d?l^D)4s4p^F 

6 

110672.3 

111988 

1.19 

112402 

1.56 

114951 

3.87 

19 

3d^(^D)4s4p*D 

8 

110867.2 

110683 

-0.17 

111004 

0.12 

113435 

2.32 

20 

3d^(^D)4s4p*D 

6 

111743.0 

111093 

-0.58 

111453 

-0.26 

113950 

1.97 

21 

3d^CD)4s4p^D 

4 

111994.3 

112044 

0.04 

112375 

0.34 

114833 

2.53 

22 

3d®(®P)4s4p®P 

8 

112409.7 

112763 

0.31 

113204 

0.71 

115809 

2.91 

23 

3d®(®P)4s4p^P 

2 

112534.9 

113409 

0.78 

113702 

1.04 

115135 

2.31 

24 

3d®^P)4s4p®P 

2 

113492.9 

112272 

-1.08 

112571 

-0.81 

116293 

2.46 

25 

3d^CD)4s4p^P 

4 

113499.2 

113479 

-0.02 

113768 

0.24 

116476 

2.62 

26 

3d^l^D)4s4p^D 

4 

114045.03 

114126 

0.07 

114440 

0.35 

117017 

2.61 

27 

3d?l^D)4s4p^D 

6 

114833.95 

114759 

-0.07 

115085 

0.22 

117659 

2.46 

MSD 




730 


1146 


3809 



Note. — NIST - experimental energies from NIST database; CItro - our HF data with a complete Cl 
expansion using TRO; CItrq - our HF data from the reduced Cl expansion calculation; CItrq - our HF data 
from the reduced Cl expansion calculation; 5 Ttro - percentage deviations of HF data from the experimental 
energies; dT^Ro " percentage deviations of the reduced Cl data from the observed energies; “ percentage 

deviations of the reduced Cl data from the observed energies. 
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gfci 


Fig. 1.— Comparison of the weighted oscillator strengths gf determined in the CItro to those 
determined in ^^tro approaches for Zn II. Dashed lines show the 25% deviation limits. 







Table 2. 


Transition line strengths S (in a.u.) for Zn II determined in the CItro approximation. 


Data 

Type 

Ni 


S 

S 

El 

1 

2 

4.199E-f00 

S 

El 

1 

3 

8.401E-f00 

S 

M2 

1 

3 

1.405E-f02 

S 

E2 

1 

4 

1.525E-f00 

S 

E2 

1 

5 

9.898E-01 

S 

Ml 

1 

6 

1.115E-06 

S 

E2 

1 

7 

6.011E-f01 

S 

E2 

1 

8 

9.013E-f01 

S 

El 

1 

9 

l.OOlE-02 

S 

El 

1 

10 

1.593E-02 

s 

M2 

1 

10 

2.201E-02 

s 

E3 

1 

11 

9.030E-06 

s 

M2 

1 

11 

1.250E-f01 

s 

El 

1 

12 

2.970E-02 

s 

M2 

1 

12 

2.061E-f00 

s 

E2 

2 

3 

1.278E-f02 

s 

Ml 

2 

3 

1.300E-f00 

s 

E3 

2 

4 

3.139E-f00 

s 

M2 

2 

4 

8.302E-03 

s 

El 

2 

5 

1.196E-02 

s 

M2 

2 

5 

7.790E-03 

s 

El 

2 

6 

2.822E-f00 

s 

El 

2 

7 

1.317E-f01 

s 

M2 

2 

7 

8.199E-f00 


Note. — The first column describes the 
transition data type (S: line strength S, 
A: transition probability A). The second 
column describes transition line type, Ni 
is for the lower level index, denotes 
the upper level index. 


Note. — (This table is available in its 
entirety in a machine-readable form in the 
online journal. A portion is shown here for 
guidance regarding its form and content.) 
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E2, E3, Ml, M2). The line strengths S are preferred as they do not depend explicitly on the 
transition energy AE. For the further use, oscillator strengths / or transition probabilities A 
can be derived using well-known relations. A complete table of transition parameters is available 
on-line. 


In Table we provide the transition wavelengths A, oscillator strengths gf, and A values for 


several observed lines. We note that for 5 of these lines, no gf or A values are listed in Morton 


(2003). 


Along with the radiative transition data, we have determined electron-impact excitation pa¬ 
rameters for these lines in the plane-wave Born approximation. The methods and codes for such 


calculation are described by Bogdanovich et al. (2014). We present just a sample of the collisional 
parameters in Table whereas the complete version of the table is available on-line. We note 
that the approximation adopted in our calculations can not produce highly accurate data for the 
electron-impact excitation process. Nevertheless, as the more accurate calculations for this level 
set of Zn II are absent, our data are an improvement over other more rough approximations, such 
as g-bar approximation. 


As already mentioned in Sect. 1.2, several elaborate studies on the Zn II ion electron-impact 


excitation were published, see Pindzola et al. (1991); Zatsarinny Sz Bandurina (1999); Sharma et 


al. (2011). Since their data are incompatible with our model due to a different atomic structure, a 


narrow energy range (as in Pindzola et al.| (1991); Zatsarinny Sz Bandurina ( 1999[ )), and a different 
data type, when only excitation cross sections are plotted instead of collision rates, we cannot 
include those data into our data set or perform a comprehensive comparison with our data. Some 
qualitative assesment suggests that our data deviate no more than 30% from the R-matrix or 
relativistic distorted-wave results at low electron energies for the optically allowed transitions and 
no more than 45% for the forbidden transitions. At high electron energies, these deviations decrease 
at least twofold. 


The ab initio calculation results for the level energies E, the radiative transition parameters 
- oscillator strengths gf, transition line strengths S, transition probabilities S, and the radiative 
lifetimes r for the Zn II are available from the database ADAMANT (http: //www. adamant. tf ai. 
vu.lt/database) being developed at Vilnius University. 

We note in passing that the theoretical calculations of the transition wavelengths can never 
be as accurate as the experimental data. We have therefore used the experimental A values in our 
Cloudy simulation runs, as has been done in previous studies. The oscillator strengths or transi¬ 
tion probabilities need to be corrected for the difference between the experimental and theoretical 
energies. We have indeed corrected the transition rates by using the experimental level energies to 
determine the radiative transition parameters (e.g., the transition probabilities A or the oscillator 
strengths gf). 
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Table 3. Transition vacuum wavelengths (in A), the upper levels, oscillator strengths gf and 
transition probabilities A (in s“^) for Zn II, including some commonly observed lines. 


A (A) 

Upper level 

gf 

A 


2062.6604 

3d^“4p "Pi/2 

6.18E-1 

4.85E+8 

2026.1370 

M^Hp^P3/2 

1.26E+0 

5.12E+8 

986.5237 

3d^%p 2Pi/2 

3.08E-3 

1.06E+7 

984.1414 

ZS%p^P3/2 

4.92E-3 

8.46E+6 

949.4630 

3d®4s4p '‘P 3/2 

9.50E-3 

1.76E+7 

938.7130 

3d®4s4p '^Pi /2 

3.65E-3 

1.38E+7 

923.9760 

ZdHsAp ‘‘P 3/2 

9.17E-6 

1.79E+4 


Note. — All these lines originate from the ground 
level 3d^°4s ^ 51/2 (<7 = 2). The statistical weights 
g — 2J + 1 of the upper levels are given in Table 







Table 4. Table 4: Effective collision strengths T for the electron-impact excitation of Zn II at 14 selected temperatures 

determined in the plane-wave Born approximation. 
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3. Cloudy calculations 
3.1. Application to DLAs and Sub-DLAs 


One of the motivations behind our interest in the Zn II atomic data stems from the somewhat 
surprising results obtained from observations of Zn II lines in DLAs and sub-DLAs. We now 
describe these results and discuss the implications of our revised atomic data for the evolution of 
DLAs/ sub-DLAs. 


Most models of chemical evolution predict a near-solar mean interstellar metallicity for galaxies 
at redshifts z ~ 0 [e.g., Pei et ah (1999); Somerville et al. ( 2001[ )]. Surprisingly, DLAs at z < 1.5 
evolve little if at all, with metallicities far below the predictions of models based on the cosmic 


star formation history [e.g., Pettini et al. () 

.999) 

Kulkarni & Fall 

(2002) 

Prochaska et al. 

(2003 

2007 

); 

Khare et al. 

(2004 

); Kulkarni et al. 

(2005 

2007 

2010 

); 

Peroux et al. 

(2008 

)]. The DLA 


global mean metallicity shows some evolution at high z, but seems to reach only about 20% of 
the solar level hy z = 0. These results appear to contradict the near-solar mass-weighted mean 
metallicity at z ~ 0 predicted by most models and observed in nearby galaxies. Equally surprising, 
a significant fraction of sub-DLAs appear to be highly metal-rich (near-solar or super-solar), even 


at redshifts z > 2 (e.g., Akerman et al. ( 

2005); Peroux et al. (2006a|b 

2008), Prochaska et al. 

(2006) 

Meiring et al. (2007 

2008 

2009a|b 

; Kulkarni et al. 

(2007 

2010 

); 

N'estor et al. 

(2008) 

Som 


et al. (2013, 2014)). The super-solar metallicities observed in many sub-DLAs (a large fraction of 


which are derived from Zn) are particularly striking, because no local counterparts to such systems 
are known among normal galaxies. 


Given these surprising results, it is natural to ask to what extent the results are affected by 
ionization of the absorbing gas. Ionization corrections are likely to be especially important for 
sub-DLAs and low-A h i DLAs. Cloudy calculations using existing atomic data suggest that the 
low-A h/ sub-DLAs can be significantly ionized, but give relatively small ionization corrections 


HM + Kurucz Spectrum HM + Kurucz Spectrum 

log N„ = 19.55, = 2.139 log N„, = 19.55, = 2.139 



Fig. 2.— Comparison of the ionization corrections and ionization parameters derived for the sub- 
DLA toward Q1039-27I9 with Cloudy version C13.02 (a) before and (b) after incorporating our 
revised atomic data for Zn 11. 




































































































































































HM + Kurucz Spectrum 

logN|jj= 19.35, = 2.058 


HM + Kurucz Spectrum 

logN|^i= 19.35, z_|^_ = 2.058 



Fig. 3.— Comparison of the ionization corrections and ionization parameters derived for the sub- 
DLA toward Q2123-0050 with Cloudy version C13.02 (a) before and (b) after incorporating our 
revised atomic data for Zn 11. 


HM + Kurucz Spectrum 


HM + Kurucz Spectmm 


log = 21.10, = 2.908 log = 21.10, = 2.908 



Fig. 4.— Comparison of the ionization corrections and ionization parameters derived for the DLA 
toward Q2342+34 with Cloudy version C13.02 (a) before and (b) after incorporating our revised 
atomic data for Zn II. 
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to the abundances (e(X) = [X/H]totai — [XII/[HI]) over the range of ionization parameters {U = 
n^lnu = ‘h 9 i 2 /cnH, where <^912 denotes the flux of radiation with \\u > 13.6 eV and uh denotes 
the total gas density) allowed by observed ratios such as A1 lll/Al II, S III/S II, or Fe IIl/Fe II [e.g., 


Dessagues-Zavadsky et al. 

(2003) 

Dessagues-Zavadsky et al. (2004); |Meiring et al. (2007, 2009b 

Som et al. 

(2013 

2014 

)]• 


To illustrate the effect of ionization, we now show the Cloudy photoionization modeling calcu¬ 
lations for a few illustrative absorbers. As typical examples of sub-DLAs, we choose the 2 ; = 2.139 
system toward Q1039-2719 and the 2 ; = 2.058 system toward Q2123-0050. Both of these sub-DLAs 
have log Nm = 19.35. As an example of DLAs, we use the 2 ; = 2.908 absorber toward Q2342-I-34 
with logAni = 21.10. We assume that the ionizing radiation incident on the absorbing cloud is 
a combination of the extragalactic UV background and an 0/B-type stellar radiation field. We 
adopt the extragalactic UV background from Haardt Madau| (1996) and Madau et al. (1999), 
evaluated at the absorber redshift. The 0/B type stellar radiation field corresponds to a Kurucz 
model stellar spectrum for 30,000 K. The incident radiation field was taken to be a mixture of the 


extragalactic and 0/B type stellar radiation fields in equal parts. Schaye (2006) has suggested that 


the contribution from local sources to the ionization of DLA absorbers may be significant when 
compared with the contribution from the extragalactic background ionizing radiation. Addition¬ 
ally, our simulations include the cosmic microwave background at the absorber redshift, and the 
cosmic ray background. However, we do not include the radiation from local shocks produced by 
supernovae, white dwarfs, or compact binary systems. 


For each absorber, we first ran grids of photoionization models using Cloudy version C13.02 
, by varying the ionization parameter from 10“® to 1. The models were made 
to match the observed H I column density and the observed metallicity based on Zn II. Constraints 
on the ionization parameter were estimated by comparing the observed values of the column den¬ 
sity ratios for various ions with the values calculated from our simulation grids. The ionization 
parameter thus estimated was used to obtain the ionization correction values to the abundances. 
In particular, we used column density ratios of adjacent ions of the same element, because they 
provide more reliable observational constraints than the ratios involving different elements, since 
the latter may be affected by differences in dust depletion or in nucleosynthesis. For the two il¬ 
lustrative absorbers discussed above toward Q1039-2719 and Q2123-0050, we estimate ionization 
parameters log U of -2.87 and -2.35, respectively. The corresponding estimates of ionization cor¬ 
rections for Zn II are 0.48 dex for the absorber toward QI039-27I9 and 0.59 dex for the absorber 
toward Q2123-0050, respectively. For S II, the corresponding ionization corrections are -0.20 dex 
and -0.30 dex, respectively. Figures and[^ show the ionization corrections for several elements 
as a function of the ionization parameter for these models, and the range of U allowed by the 
observed ion ratios in these absorbers. For the DLA toward Q2342-I-34, we estimate ionization 
parameter log 7/ of -3.41 using Al III/SIII and Zn II ionization correction of -1-0.16 dex (Fig. [^). 

To assess the effect of our atomic data calculations, we next ran the revised Cloudy models that 
include our improved Zn atomic data for the illustrative DLA and sub-DLA absorbers discussed 


Ferland et al.| (2013)] 
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above. Figures and[^ show the results for the ionization correction as a function of ionization 
parameter, and the range of U allowed by the observed ion ratios. The ionization parameters are 
-2.90 and -2.38, respectively, for the absorbers toward Q1039-2719 and Q2123-0050. The corre¬ 
sponding ionization corrections for Zn II are 0.45 dex and 0.54 dex, respectively, i.e. lower than 
those obtained with Cloudy version C13.02 by 0.03 dex and 0.05 dex, respectively. For the DLA 
toward Q2342-I-34, we estimate ionization parameter log?7 of -3.40 using A1 Ill/Si II and Zn II 
ionization correction of -1-0.16 dex (Fig. |^). 

Combining the difference in ionization correction with the difference in A^znii implied by the 
revised oscillator strengths (which are lower by 0.10 dex), the revised [Zn/H] values for the above 
sub-DLA absorbers would be lower by 0.13—0.15 dex compared to the values based on the previously 
available atomic data. For the DLA absorber, the revised [Zn/H] value would be lower by 0.10 dex 
compared to the values based on the previous atomic data. 

We note, however, that the estimates of ionization corrections are sensitive to the adopted 
values of dielectronic recombination rates, which are unknown for ions of most elements in the 
third row of the periodic table and beyond, including key elements such as Al, Fe, Zn. In future 
papers, we plan to address the dielectronic recombination rates. 


3.2. Emission lines 

The atomic data described in this paper will be part of the next release of Cloudy. For 
reference, Figureshows the lowest energy levels and indicates some of the stronger lines with their 
air wavelengths given in A. We know of no calculations of the emission spectra that demonstrates 
its diagnostic power, so we do representative calculations here. 

Figure shows Zn II emission spectra at a temperature of 10^ K and electron densities of 1 
and 10^*^ cm“^. The gas was assumed to be composed entirely of Zn"*" with no incident SED, so 
the emission is entirely due to electron impact excitation. The strongest lines, as expected, are 
the 3d^°4p ^ 73 / 2 , 1/2 ~ 3d^‘^4s doublet at AA 2025.48A, 2062. OOA. The next stronger lines are 
3d®4s^ ^ 7 ) 5/2 — 3(i^‘^4p ‘^Pz/2 7478.82A and 3(i®4s^ “^^3/2 ~ 37^*^4p ‘^Pi/2 at 5894.37A. These 

lines are considerably fainter but are important because they can be detected with large ground- 
based instruments. The same is true for the El doublet lines 3(i^‘^5s ‘^Si /2 — 3d^^4:p ^7’i/2,3/2 at AA 
2501. 99 A, 2559. 95 A and the doublet of the lines at AA 2064.23A, 2099.94A, representing the El 
transitions — Sd^^ip. 

Figure shows several possible Zn II temperature indicators. To do this, the kinetic temper¬ 
ature of a pure Zn'*' gas with an electron density of I cm“^ was varied over a wide range and the 
emission ratios suggested by Figure [^plotted. 

The upper panel shows the ratio of the two optical lines relative to the sum of the two strongest 
UV lines. This is an indicator with a wide dynamic range although the plot also shows that the 
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Fig. 5.— The diagram of the lowest levels for Zn II, showing the even-parity configuration 

levels and the odd-parity configuration ‘Ad^AsAp levels. The strongest 

lines originating from the “Ad^^Ap — 3(i^'^4s transition (solid lines) and from the A(^As^ — Ad}^Ap 
transition (dashed lines) are also plotted. Other strong doublet El lines, representing transitions 
3d^°5s — Ad}^Ap (dotted line) and Ad}^Ad — Ad^^Ap (dash-dotted line) are also shown. 
















log intensity [arbitrary units] 
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Wavelength [gm] 


Fig. 6.— The Zn II strongest line and continuum emission. The gas has only Zn"*" at 10^ K and 
is computed at two densities. The vertical scale is adjusted so that the two spectra lie near one 
another, so the vertical scale is arbitrary. Both Zn II recombination and brews emission are present. 
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Fig. 7.— Several Zn II temperature indicators are shown. 
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optical lines are far fainter than the UV transitions. The lower panel shows the ratio of the two 
optical lines. These have the advantage of being detectable by ground-based instrumentation. 


4. Conclusions 


Our estimates of the oscillator strengths for the key Zn II absorption lines at 2026.14, 2062.66 
A are higher than the previous values by 0.1 dex, implying the Zn II column densities inferred 
from these lines to be lower by 0.1 dex. Moreover, the sub-DLA ionization corrections for Zn II 
would be lower by ~ 0.03 — 0.05 dex, as discussed in section 3.1 Thus, the logarithmic metallicities 
inferred from these lines would be lower by 0.13 — 0.15 dex for sub-DLAs and by 0.1 dex for DLAs, 
compared to past studies. While differences of this amount are significant, they are not adequate 
to explain why the sub-DLA metallicities are so much higher than that of DLAs. 


Using Cloudy simulations, we have demonstrated some astrophysical applications of our atomic 
data calculations, and the predictions for many emission and absorption lines. One can compare 
such predictions with the observed line strengths to obtain improved constraints on the chemical 
composition and physical properties of the Galactic and extragalactic ISM. Past observations of 
Zn II absorption in the Galactic ISM and the DLAs/sub-DLAs have targeted the AA2026.I4, 2062.66 
lines. Our new calculations show that these are indeed the strongest observable lines in the com¬ 
monly accessible wavelength region. Our complete set of Zn II oscillator strengths (available online) 
contains many UV transitions. Although most of these transitions are weak, some of them may be 
detectable in high-S/N spectra obtained with the extremely large telescopes of the future. 
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